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Abstract 

Anomalous isotope shifts in the chain of light Hg isotopes are investigated by using the Hartree-Fock- 
Bogoliubov method with the Skyrme SIII, SkI3 and SLy4 forces. The sharp increase in the mean-square 
radius of the odd mass 181_185 Hg isotopes is well explained in terms of the transition from an oblate to a 
prolate shape in the ground state of these isotopes. We discuss the polarization energy of time-odd mean-field 
terms in relation to the blocked level by the odd neutron. 
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1 Introduction 

The sharp increase in the mean-square charge radius is observed in the odd mass 181 - 185 Hg isotopes with respect 
to their even neighbors [1]. A hint that a very different configuration may constitute the ground state of these 
nuclei was given by an isomer shift measurement of 185 Hg [2], where they found an excited state of 185 Hg which 
has a charge radius consistent with the charge radius of heavier neighbors. The sharp increase in the charge 
radius was interpreted as the transition from a nearly spherical to a deformed shape in the ground state of the 
odd mass 181 - 185 Hg isotopes. 
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Frauendorf and Pashkevich [3] calculated the deformation energy of Hg isotopes by means of Strutinsky's 
shell correction method. They found that the light Hg isotope has an oblate and prolate minimum of nearly the 
same depth at the respective deformations e = —0.12 and 0.22. They discussed a mechanism which causes the 



that the neutron level density is very high in the almost spherical oblate minimum, whereas the neutron level 
density is low in the well deformed prolate minimum. They concluded that the loss of pairing energy due to 
the blocking of the Fermi level by the odd neutron is larger in the oblate minimum than in the prolate one, 
which gives rise to the oblate-prolate transition in the ground state of these isotopes. As the prolate minimum 
corresponds to deformation twice as large as the oblate one, the sharp increase in the mean-square radius is 
observed in the odd-mass light Hg isotopes. 

There are several mean field calculations of Hg isotopes. However, all of them are limited to even isotopes 
[4, 5, 6, 7]. In order to see if the mechanism of Frauendorf and Pashkevich works well, we perform Hartree- 
Fock-Bogoliubov (HFB) calculations of Hg isotopes with effective interactions SIII [8], SkI3 [9] and SLy4 [10]. 
We calculate Lipkin-Nogami corrections [11, 12] to both mean field and pairing energies [13, 14]. We investigate 
the influence of the time-odd mean-field terms as well as the blocking of a level by the odd neutron on the 
oblate-prolate shape transition. 

In Section 2, we present the Skyrme HFB method with Lipkin-Nogami corrections to describe both even 
and odd Hg isotopes. In Section 3, we present calculations with effective interactions SIII, SkI3 and SLy4. 
Conclusions are given in Section 4. 

2 Method of calculation 

2.1 Skyrme Hartree-Fock-Bogoliubov method 

The HFB equation is derived by the variation of energy with respect to the quasi-particle vacuum [15], 
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Here the quasi-particle vacuum |<j>) is given by 
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for all quasi-particle states k. In the present study, we employ Skyrme forces in the particle-hole channel and 
the zero-range density-dependent pairing force in the particle-particle channel, 



The HFB energy is then given by 

where H is the Skyrme energy density [16, 17] and A(r) is the gap field, 



(5) 



(6) 



A W = |(l-^]M,, 



(7) 



The pairing strength Vq is determined from the odd-even mass difference and the density parameter p c is taken as 
the saturation density of symmetric nuclear matter which is obtained by using the respective Skyrme parameters. 
The particle density p(r) and the anomalous density n{r) are written as 



p(ra, ru) = ^ p a p<p a (r, cr)(p* (r, a'), 
n(ra, ra') = ^ K al3 ip a (r, a)<pp(r, a') 

afi 

with the density matrix and pairing tensor, 

k 

K a p = (^\c p c a \^) = V* k U pk . 
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Variation of Eq. (6) with respect to p a/ 3 and n a p yields the HFB equation, 
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Diagonalization is carried out in a deformed oscillator basis with axial symmetry [18]. A Greek index in Eq. (12) 
denotes a set of quantum numbers {n r n z AS} of a deformed oscillator wave function. 



In the present study we will restrict our basis space as 

E a = (2n r + |A| + l)Ku)j_ + [n z + hu z < ^N + hu (13) 

with N = 18. The oscillator constant is held fixed as Tiu = h{J\uj z ) 1 l ? ' = 1.21 x ^A^ 1 / 3 . The deformation 
parameter q = uj±/uj z is varied to obtain the maximum deformation energy in the oblate and prolate minimum, 
respectively. The chemical potential A is adjusted by the particle number condition, (N) — N. Using the 
occupation amplitude N k = J2 a VakYak and quasi-particle energy E k , we calculate an auxiliary single-particle 
energy and pairing gap in order to find the optimum value of A [19], 

e k = -E k (2N k - 1) + X, (14) 



A fc = 2E ky /N k (l~N k ). (15) 

As in Ref. [20, 21], the cut-off procedure of quasi-particle states is imposed on the auxiliary single-particle 
energy. The value of e max = 60 MeV is assumed in the present study. 

2.2 Description of an odd nucleus 

The above formalism can be straightforwardly extended to the case of one quasi-particle state in odd mass nuclei. 
All the equation retain their form. Only the density matrix p a p and pairing tensor n a p should be modified due 
to the blocked quasi-particle /3^|$) in the following way [22, 23]; 

Pa/3 = p a p + U ak Up k - V* k V 0k , (16) 

a a /3 = Had + u ak Vp k - v* k Up k - (17) 

We introduce a reference state |$hfbe) and its energy expectation value E HFBE . The state is constructed as 
an even vacuum without quasi-particle creation and without breaking time-reversal invariance but with an odd 
average particle number. The difference of energy £? HFB — £?hfbe j g ^.^g com p 0sec i f the loss of pairing energy 
due to the blocking of the Fermi level and the polarization energy of the mean field by the presence of an odd 
neutron [24]. 

A proper description of an odd nucleus by mean field theories requires to break the time-reversal symmetry 
[17]. When the time-reversal symmetry is broken, the Skyrme energy density H is written by the sum of kinetic 
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energy and potential energy of isoscalar 7i and isovector 7ii terms; 

h 2 

Tt = - — tq+Ho+Hi, 
2m 

H t = C p tP 2 + C s t sl + Cf p p t & Pt + Cf s SfAs t 

+C T t { Pt T t - j 2 ) + C/(J t 2 - a t ■ T t ) + C? J [p t V ■ J t + s t ■ (V x jt)]- (18) 

The coupling constants C are expressed by the Skyrme force parameters [25]. The isospin index t can have 
values or 1. The isoscalar and isovector parts of the particle densities p t are defined as 

P 0=Pn+Pp, Pl=Pn-Pp, (19) 

and analogous expressions are used to define other densities. It should be noted that since we consider the 
Skyrme force as the density-dependent two-body interaction, the time-odd terms C t s in Eq. (18) are different 
from the original expressions [17] which were derived from the three-body interaction. 

The characteristic combinations of (ptn — j 2 ), (J 2 — s t ■ T t ) and [p t V • J t + s t • (V x j t )] in Eq. (18) are 
the results of the local gauge invariance of the Skyrme force [25] . The time-odd coupling constants are uniquely 
determined from the original Skyrme-force parameters and therefore the time-odd mean field has no adjustable 
parameters once the time-even mean field is given. However, as discussed by Dobaczewski and Dudck [25], one 
may consider the energy density to be a more fundamental construction than the Skyrme interaction itself. In 
such a case, all 20 coupling constants in Eq. (18) can be treated and adjusted independently. However, the 
gauge-invariant conditions restrict the values of six time-odd coupling constants and leave the freedom to modify 
the values of C t s and C As . We are going to vary C t s and C As to see the influence of these terms on the oblate- 
prolate transition in the odd-mass 181 ~ 185 Hg isotopes. We will omit the terms (J 2 — s t ■ T t ) because J\ are 
omitted in the parameterization of the present Skyrme forces. 

2.3 Lipkin-Nogami corrections 

In light Hg isotopes, the neutron level density is high in the oblate minimum, while it is low in the prolate 
minimum. Because of the difference of level densities in the respective minima, it happens that neutrons in the 
oblate minimum are superconductive, whereas they are normal in the prolate minimum within the same isotope. 
The situation is extremely inconvenient because we want discuss the energy difference between the oblate and 
prolate minimum. Besides, a sudden phase-transition from the normal to superconductive state or vice versa 
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at certain neutron numbers may be unphysical and it may be due to the defect of the number non-conserving 
pairing scheme. In the present study we employ the approximate number-projection method of Lipkin and 
Nogami (LN) [11, 12] and calculate its corrections to both mean field and pairing energies [13]. The number 
projected energy is expressed as 

£hfbln = £hfb + Ai(JV - (AO) + A 2 (iV 2 - (A> 2 )) + • • • , (20) 

where Ai and A2 are constant parameters [14]. The second order expansion is considered in the Lipkin-Nogami 
method. Variation with respect to the density matrix and pairing tensor as in the case of the HFB equation yields 
the HFBLN equation applicable to both even and odd nuclei (hereafter simply called HFB). The parameter Ai 
is given by the particle number condition, whereas the parameter A 2 is held fixed during the variation and is 
determined after the variation from the additional condition, 

|0=exp(i£A>2 )|$>, (21) 
dl($\H\0 



A 2 = 



(22) 



£=0 

Equation (22) is calculated by the following density matrix and pairing tensor [13]; 



Pal = {^ f3 c a \C) = (fi + 2ii{l-p)p) afil (23) 



JO 

K af3 



<$|c j8 c a |£> = (K + 2i£(l-p)/ S ) ai9 , (24) 
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a/3 — \^l°Q°/3 



^\cicl\0 = (K*-2i^*pU. (25) 



For an odd nucleus, the density matrix and pairing tensor of Eqs. (16) and (17) are used on the right hand side 
of the above equations [14]. 



3 Results and discussion 

We have analysed the oblate-prolate transition in the odd-mass light Hg isotopes using the Skyrme Hartree- 
Fock-Bogoliubov method with effective interactions SHI, SkI3 and SLy4. Shape transition occurs only in the 
odd-mass isotopes as the result of a subtle balance between the time-even mean-field energy, the loss of pairing 
energy by the odd neutron and the time-odd mean-field energy. We have examined effective interactions SIII, 
SGII [26], SkM* [27], SKX [28], Skll, SkI2, SkI3, SkI4 and SkI5 [9] and SLy4 whether they can predict the 
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oblate-prolate shape coexistence in light Hg isotopes. We have found that both SIII and SkI3 interactions have 
desirable isospin dependence and deformation energies. 

We have assumed the zero-range density-dependent pairing interaction. The pairing strength Vq was adjusted 
to reproduce the one-particle separation energy [29] . The density parameter p c is taken as the saturation density 
of symmetric nuclear matter which is obtained from Skyrme parameters. Adjusted values are given in the 
respective figure captions. 

We have investigated the effect of time-odd mean-field terms on the HFB energy, first by using the time-odd 
coupling constants derived from the Skyrme force parameters and second by switching on and off the coupling 
constants C t As and Cf, respectively. These coupling constants in C-representation [25] are given in Table 1. 
They vary very much with different Skyrme forces, e.g., the coupling constant Cq(p = 0) is 14, 84 and -208 for 
the respective SIII, SkI3 and SLy4 forces. 

The polarization energy depends strongly on the blocked level by an odd neutron. Even the sign of energy 
changes with different choice of the blocked level. In the present study we simply choose as the blocked level the 
lowest quasi-particle state that is obtained from the HFBE calculation. 

Mercury isotopes are very interesting because the shape of the ground state varies very much as the neutron 
number decreases from the closed shell. Figure 1 shows the deformation energy of even Hg isotopes calculated by 
the constraint SIII-HF method. We see the spherical shape in the neutron-closed 206 Hg. The oblate minimum 
with deformation (3 2 ~ —0.17 then develops in 200 Hg. A rather complicated prolate minimum appears in 190 Hg. 
Finally the deep prolate minimum with deformation 02 ~ 0.29 develops in 184 Hg. The oblate-prolate shape 
coexistence in light Hg isotopes has been predicted by many authors, e.g., [3, 4, 5, 6, 7, 30, 31]. 

Figure 2 shows the SIII calculation of the oblate-prolate energy difference. The HF calculation (white circle) 
predicts prolate deformation in the ground state of all Hg isotopes. The HFBE calculation (black circle), the 
fully paired state without the blocking effect, predicts oblate deformation in the ground state of all Hg isotopes. 
The shape changes due to the large gain of pairing energy in the oblate minimum because of the high level 
density as compared to the level density in the prolate minimum. In the SIII- HFB calculation (white square), 
on the other hand, the shape of odd mass 181 ~ 185 Hg isotopes returns to the prolate shape. This is because the 
loss of pairing energy due to the blocking by the odd neutron is also large in the oblate minimum than in the 
prolate one. The effect of the time-odd mean fields on the oblate-prolate transition is weak and a few tenths of 
the blocking effect. 
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Figure 3 shows the SIII-HFB calculation of (a) the mean-square charge radius and (b) one-neutron separation 
energy of 18 °- 206 Hg isotopes. The sharp increase in the mean-square radius of the odd-mass 181 - 185 Hg isotopes 
is very well explained by the calculation. The one-neutron separation energy is also well reproduced in the 
calculation. This is so because the depth of the pairing force has been adjusted to reproduce the observed 
separation energies. It is noted that the white circle in Fig. 3(b) is a derived value from systematic trends [29] 
and not an experiment. 

As is well known, a proper description of odd nuclei by mean field theory requires to break the time-reversal 
symmetry and introduces the time-odd components (TOC) in the mean field. The effect of an odd neutron on 
the HFB energy is therefore twofold; (a) The loss of pairing energy due to the blocking of the Fermi level by the 
odd neutron, and (b) the time-odd mean-field energy arising from the blocked level. The time-odd mean fields 
have been studied both in the analysis of super-deformed rotational bands [25] and the odd-even nuclear mass 
difference [24]. 

In order to see the blocking energy (the loss of pairing energy) and the polarization energy (the time-odd 
mean-field energy) separately, we define the blocking energy £; block as the energy difference between E HFB with 
no TOC and £: HFBE . We also define the polarization energy E po1 as the energy difference between E UFB and 
£ HFB with no TOC, i.e., 

E hlock = £ HFB (no TOC) -E HFBE , (26) 

£ P ol = E HFB _£HFB (no T0C) (27) 

Figure 4 shows the SIII-HFB calculation of (a) the blocking energy in the oblate and prolate minimum, (b) 
the polarization energy in the oblate minimum, (c) the polarization energy in the prolate minimum and (d) the 
oblate-prolate difference of HFB energy. In Fig. 4(a), the blocking energy in the oblate minimum is larger by 
^200 keV than the one in the prolate minimum. This results from higher neutron level-density in the former 
than in the latter. The blocking energies are reversed in 189 Hg because of the small deformation (3 = 0.08 in the 
prolate minimum of this isotope and the level density in the minimum is high accordingly. 

As mentioned in Section 2.2, the C t s and C t As terms are free from the time-even parameterization. To 
investigate their effects on the HFB energy, we have further divided the polarization energy into contributions 
with (i) the full time-odd terms, (ii) the time-odd terms without Cf, (hi) the time-odd terms without C t As , and 
(iv) the time-odd terms without both C t s and C t As . In the following we assume simultaneous modifications of 
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the isoscalar and isovector coupling constants of a given species. 

Figures 4 (b) and (c) display the contributions to the polarization energy. The polarization energy with full 
time-odd terms (white circle) is quite small. It is attractive and less than -50 keV in the oblate minimum, while 
it is repulsive and less than 50 keV in the prolate minimum. When the Cf terms are omitted (black circle), 
the polarization energy becomes attractive and more than -250 keV in the oblate minimum and it is attractive 
around -150 keV in the prolate minimum. The C t s terms have repulsive contribution to the HFB energy When 
the C As terms are omitted (white square), on the other hand, the polarization energy becomes repulsive and 
^50 keV in both minima. The C As terms have large attractive contributions to the HFB energy. When both 
Cf and C t As terms are omitted (black square), the polarization energy becomes small. The effect of the C t s and 
C t As terms is opposite and cancels with each other. 

Figure 4(d) shows the difference of energy in the oblate and prolate minimum. The figure shows the calculation 
with the full time-odd terms (white circle), without the C t s terms (black circle), without the C As terms (white 
square), without both C t s and C t As terms (black square) and without all time-odd terms (white triangle). When 
both CI and C t As terms are omitted, the oblate shape is predicted in 181 Hg, while the prolate shape is calculated 
in 185 Hg with marginal energy difference between the oblate and prolate minima. 

Table 2 summarizes the oblate-prolate difference of the HFBE energy A£? HFBE , the blocking energy AE hlock , 
the polarization energy AE po1 and the sum of these contributions A£? HFB . The polarization energy AE po1 is 
calculated with four different time-odd mean fields. They are the HFB calculation (i) with full time-odd terms, 
(ii) without the C t s terms, (iii) without the C t As terms and (iv) without both C t s and C t As terms. 

The energy difference A_E HFBE shows that the minimum is found at the neutron number TV = 102 and 
the oblate minimum gets deeper than the prolate one as the neutron number is away from N = 102 in both 
directions. For the blocking energy E hlock , it is repulsive and larger by 400 keV in the oblate minimum than in 
the prolate one. The blocked level is labeled with the quantum number 0, the component of neutron angular 
momentum along the symmetry axis. 

For the odd mass 181 - 185 Hg, the blocked level has fl = 1/2 in the oblate minimum, while it has il — 5/2 
or 7/2 in the prolate minimum. In the calculation with full time-odd terms (case (i)), the polarization energy 
is attractive in the oblate minimum, while it is repulsive in the prolate minimum, leading to the the attractive 
energy difference of A_E po1 . The attractive Ai? po1 weakens the repulsive Ai? block and prevents the oblate-prolate 
shape transition in these isotopes. In the calculation without the C t s terms (case (ii)), the polarization energy 
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becomes attractive because of the absence of the repulsive C t s terms. In this case, the polarization energy 
A_E po1 almost cancels the blocking energy AE Uock . In the calculation without the C t As terms (case (hi)), the 
polarization energy is weakly repulsive because of the absence of strongly attractive C t As terms. Finally in the 
calculation without both C t s and C t As terms, the polarization energy is weakly attractive in both minima. The 
effects of the C t s and C t As terms on the polarization energy will be discussed more generally in the last paragraph 
of this section. 

Figure 5 shows the results of the SkI3-HFB calculation. We had the difficulty of convergence in this calculation 
when the full time-odd components are considered. The difficulty arises mainly due to the strong coupling- 
constant C As derived from the SkI3 parameters (see Table 1) and partly due to the large HF basis of the present 
study (N = 18). The difficulty occurs when the energy splitting between the time- reversed partners becomes 
large and a few of the quasi-particle states have negative energy. Since we have no way to correct this difficulty, 
we simply omit the C As terms and display the results in Fig. 5. The blocking energy E hlock in Fig. 5(a) is similar 
to the one in the SIII-HFB calculation. Figures 5(b) and 5(c) show the polarization energy (iii) without C t As 
terms (white square) and (iv) without both C As and Cf terms (black square). The Cf terms have repulsive 
contributions of ~T50 keV in the oblate minimum while it is ~200 keV in the prolate minimum. In Fig. 5(d), 
the oblate-prolate transition is very well explained when the terms C As are omitted. It is also well explained by 
omitting the whole time-odd components. 

Table 3 summarizes the results of the SkI3-HFB calculation. The energy difference AE KFBE is nearly the 
same for even and odd 180 ~ 185 Hg isotopes. This fact helps very much in reproducing the oblate-prolate transition 
in the odd mass 181 ~ 185 Hg isotopes. The SkI3 force is very well suited to describe the mean fields of light Hg 
isotopes. 

Figure 6 shows the results of the SLy4-HFB calculation. The blocking energy £; block in the odd-mass 
181-185 Hg is larger by ~200 keV in the oblate minimum than in the prolate one. However, their magnitudes are 
reversed at the neutron number N = 107 and N = 109. The reason is the same as the SIII calculation. The 
equilibrium deformation in the prolate minimum is reduced to j3 = 0.12 in these isotopes. Figures 6(b) and 6(c) 
show the terms have repulsive contributions to the polarization energy. In addition, the terms C% and C As 
have opposite signs and nearly cancel with each other. 

The oblate-prolate energy difference in Fig. 6(d) shows that the SLy4 force is not very well suited to describe 
deformation properties of light Hg isotopes. Here we also had the difficulty of convergence when the terms are 
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not included in the calculation. The difficulty may be due to the C As terms. The repulsive C t s terms cancel the 
attractive C As terms. Without the C t s terms, the attractive C t As terms are quite strong and cause the difficulty 
of convergence. 

Table 4 summarizes the results of the SLy4-HFB calculation. The energy difference Ai? HFBE shows that 
the SLy4 force favors the oblate shape in light Hg isotopes much more than the SIII and SkI3 forces. This fact 
results in the oblate-prolate transition only in the ground state of 181 Hg isotope. The oblate shape is predicted 
for all heavier isotopes. 

To see the effects of the time-odd C t s and C t As terms more closely, we have analysed the odd-mass Sn, Ba, Yb, 
Hg and Pb isotopes by the Skl3-HFB method with the basis space of No = 10. A small basis space was chosen 
in order to avoid the difficulty of convergence. Results of the SkI3-HF calculation are shown in Fig. 7. The 
polarization energy is plotted against 2 x of the blocked level. Polarization energies are calculated (a) with 
the full time-odd terms, (b) without C t As terms and (c) without both C t As and Cf terms. We may conclude that 
the C As terms are attractive and have a large effect on the blocked level with low SI, in particular £1* = l/2 ± . 
The terms CI are repulsive for levels with any value of il. The remaining time-odd terms are attractive for levels 
with high ft and are negligible for levels with low ft. The same features of the time-odd terms are observed also 
for the SIII and SLy4 interactions. For the odd-mass 181 ~ 185 Hg in Tables 2, 3 and 4, we have seen that the 
time-odd terms have attractive contributions to A_E HFB in the cases (i) and (ii), while the terms have repulsive 
contributions in the case (iv). 

4 Conclusions 

The shape and the mean-square radius of Hg isotopes have been investigated by using the Skyrme Hartree-Fock- 
Bogoliubov method. We have first examined effective interactions SIII, SGII, SkM*, SKX, Skll, SkI2, SkI3, 
SkI4, SkI5 and SLy4 to see if they can predict the oblate-prolate shape coexistence in light Hg isotopes. We have 
then examined the loss of pairing energy due to the odd neutron by employing the zero-range density-dependent 
pairing force. Lipkin-Nogami corrections to both even and odd nuclei are calculated within the HFB formalism. 
We have also investigated the effect of the time-odd C t As and C t s terms on the HFB energy in relation to the 
blocked level by the odd neutron. 

We have obtained the following conclusions. (1) The effective interactions SIII and SkI3 have desirable 
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deformation properties in predicting the oblate-prolate shape coexistence in light Hg isotopes. (2) The blocking 
energy of the pairing is the most important factor to change the shape of the odd-mass 181 ~ 185 Hg isotopes. (3) 
The effect of the time-odd terms on the prolate-oblate transition is weak and cancels the effect of the blocking 
energy by a few tenths. 

We have also investigated the effects of the time-odd terms C t s and C t As by blocking several levels near the 
Fermi surface of Sn, Ba, Yb, Hg and Pb isotopes. We have found that the terms C t As are attractive for the 
blocked levels with low 0, in particular SV = l/2 ± . On the other hand, the terms C t s are repulsive for the 
blocked levels with all values of Q. The remaining time-odd terms are attractive for the blocked levels with high 
fl and negligible for the blocked levels with low f2. 

We had the difficulty of convergence in the SkI3-HFB calculation when the C t terms are included. The 
difficulty also arises from the SLy4-HFB calculation when the C t s terms are omitted. These difficulties are mainly 
due to the strong C t As terms and partly due to the large basis space of the calculation (No = 18). For large C t As , 
the energy splitting between the time-reversed partners becomes large and a few of the quasi-particle states have 
negative energy. 
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Table 1: The time-odd coupling constants C t s and C t As for the Skyrme SIII, SkI3 and SLy4 forces. The coupling 
constants C t s are given in MeV-fm 3 , while C t As are given in MeV-fm 5 . 

Force C$(p = 0) Cg(p = p NM ) C((p = 0) Cf(p = p NM ) C As Cf s 

SIII 14.109 56.401 141.094 98.802 17.031 17.031 

SkI3 84.486 253.799 220.360 113.551 92.235 22.777 

SLy4 -207.824 153.382 311.114 99.635 47.057 14.282 



15 



Table 2: Summary of the SIII-HFB calculation in keV. Notations; fl* hl and Qp Ioi denote the blocked level 
in the oblate and prolate minimum, respectively. A£ HFBE = E™ BE - E^ BE , A£ block = E ] ^° ck - E™°f, 
AE po1 = E^°l - E^f oV and A£ HFB = £ B b FB - E™ B . Polarization contributions arc further divided into (i) the 
calculation with full time-odd terms, (ii) the calculation without Cf terms, (iii) the calculation without C t As 
terms, and (iv) the calculation without both C 4 S and C t As terms. 





N 


100 


101 


102 


103 


104 


105 


106 


107 


108 


109 






- 167 


-110 


- 68 


- 76 


- 152 


-264 


-470 


-807 


-1323 


-962 




s2 obl 


- 


1/2" 




1/2" 


- 


1/2" 


- 


1/2" 


- 


5/2- 




^prol 


- 


7/2- 




7/2- 




5/2- 




9/2+ 




9/2- 




TT'block 
^obl 


- 


1220 




1263 




1227 




1317 


- 


1186 




*W 


- 


934 




747 


- 


854 


- 


904 


- 


1412 




A -gblock 


- 


286 




516 


- 


373 


- 


413 


- 


-226 


(i) 






- 34 




- 39 


- 


- 10 


- 


- 10 


- 


7 






- 


32 




31 


- 


36 


- 


47 


- 


36 




A£p o1 


- 


- 65 




- 70 




- 46 




- 58 




- 30 




AE HFB 


- 167 


110 


- 68 


370 


- 152 


63 


-470 


-452 


-1323 


-1217 


(ii) 


^obl 


- 


-317 


- 


-328 


- 


-246 


- 


-250 


- 


-173 




<, 


- 


-141 




-102 




-142 


- 


-111 


- 


-147 




A£p o1 


- 


-176 




-226 


- 


-104 
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-139 


- 


- 26 






-167 


- 1 


- 68 


214 


-152 
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-470 


-533 


-1323 


-1214 
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ppol 
^obl 




32 




58 




79 




76 




49 




E pol 
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58 




73 




73 




61 
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- 28 









6 




3 




- 11 




AE HFB 


- 167 


147 


- 68 


440 


- 152 


115 


- 470 


-391 


-1323 


1199 
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- 32 
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- 16 




- 43 
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- 66 




- 64 




- 41 




- 41 




- 80 
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35 




32 




22 




25 




36 




AE HFB 


- 167 


210 


- 68 


472 


- 152 


131 


- 470 


-370 


-1323 


-1152 
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Table 3: Summary of the SkI3-HFB calculation in keV. Notations are the same as those in Table 2. 





N 


100 
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102 


103 


104 
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107 


108 


109 
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- 90 
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440 


-268 
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- 26 
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-147 
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73 




74 




67 




AE HFB 


- 54 
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- 93 
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- 76 
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-268 


-343 


-1159 
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Table 4: Summary of the SLy4-HFB calculation in keV. Notations are the same as those in Table 2. 
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101 
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106 
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108 


109 
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- 219 


- 276 


- 428 
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-1530 


-1601 
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- 81 
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- 42 
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- 
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Figure 1: Deformation energy of Hg isotopes calculated by using the constraint SIII-HF method. 
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Figure 2: The oblate-prolate energy difference as a function of the neutron number N of Hg isotopes. 
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Figure 3: (a) The mean-square charge radius of 180 - 206 fjg w ith respect to 198 Hg. (b) The one-neutron separa- 
tion energy as a function of the neutron number N of Hg isotopes. Experiments are taken from [1] and [29], 
respectively. White circles in (b) denote the separation energy derived from systematic trends [29]. 
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Figure 4: Results of the SIII-HFB calculation with the pairing parameters of Vq = —567 MeV and p c = 0.145 
fm~ 3 . (a) The loss of pairing energy due to the self-consistent blocking of the Fermi level by the odd neutron, 
(b) The polarization energy in the oblate minimum, (c) The polarization energy in the prolate minimum, (d) 
The oblate-prolate energy difference. Case (i): the calculation with full time-odd terms (white circle). Case (ii): 
the calculation without the C t s terms (black circle). Case (hi): the calculation without the C As terms (white 
square). Case (iv): the calculation without both Cf and C t As terms (black square). The calculation with no 
time-odd components is shown in white triangle. 
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Figure 5: Results of the SkI3-HFB calculation with the pairing parameters of Vq = —560 MeV and p c — 0.158 
fm~ 3 . (a) The loss of pairing energy due to the self-consistent blocking of the Fermi level by the odd neutron, (b) 
The polarization energy in the oblate minimum, (c) The polarization energy in the prolate minimum, (d) The 
oblate-prolate energy difference. Case (hi): the calculation without the C t As terms (white square). Case (iv): 
the calculation without both Cf and C 4 As terms (black square). The calculation with no time-odd components 
is also shown in white triangle. 
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Figure 6: Results of the SLy4-HFB calculation with the pairing parameters of Vq = —500 MeV and p c = 0.160 
fm~ 3 . (a) The loss of pairing energy due to the self-consistent blocking of the Fermi level by the odd neutron, 
(b) The polarization energy in the oblate minimum, (c) The polarization energy in the prolate minimum, (d) 
The oblate-prolate energy difference. Case (i): the calculation with full time-odd terms (white circle). Case (hi): 
the calculation without the C t As terms (white square). Case (iv): the calculation without both C t s and C 1 
terms (black square). The calculation with no time-odd components is also shown in white triangle. 
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Figure 7: The time-odd mean- field energy calculated by the SkI3-HFB method with the basis space of N = 10. 
The energy is plotted against 20 ± , twice of the component of neutron angular momentum along the symmetry 
axis, (a) The calculation with full time-odd terms, (b) the calculation without the C t As terms, and (c) the 
calculation without both C^ s and Cf terms. 
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